/ Archive Magazine CD 1995 / Archive Magazine CD 1995.iso / discs / prog_disc / volume_4 / issue_07 / help / hawk / !HawkV10 / FFt

< prev

next >

MacOS 8.1

Win98

DOS

view JSON data

view as text

+--------+-------------------------+-------------------------+--------+--------+
|00000000| 0d 00 0a 15 f4 20 3e 24 | 2e 21 48 41 57 4b 56 31 |..... >$|.!HAWKV1|
|00000010| 30 2e 46 46 74 0d 00 14 | 1c ee 20 85 20 f1 20 f6 |0.FFt...|.. . . .|
|00000020| 24 3b 22 20 61 74 20 6c | 69 6e 65 20 22 3b 9e 3a |$;" at l|ine ";.:|
|00000030| e0 0d 00 1e 07 eb 20 38 | 0d 00 28 0b 40 25 3d 26 |...... 8|..(.@%=&|
|00000040| 36 30 41 0d 00 32 0a 2a | 43 4c 4f 53 45 0d 00 3c |60A..2.*|CLOSE..<|
|00000050| 37 f4 20 46 46 74 2d 50 | 72 6f 67 72 61 6d 6d 20 |7. FFt-P|rogramm |
|00000060| 66 75 65 72 20 32 2d 64 | 69 6d 20 54 72 61 6e 73 |fuer 2-d|im Trans|
|00000070| 66 6f 72 6d 61 74 69 6f | 6e 20 6d 69 74 20 46 69 |formatio|n mit Fi|
|00000080| 6c 74 65 72 0d 00 46 0d | f4 20 31 38 2d 34 2d 39 |lter..F.|. 18-4-9|
|00000090| 30 0d 00 50 09 6e 3d 32 | 35 36 0d 00 5a 0b 6e 30 |0..P.n=2|56..Z.n0|
|000000a0| 25 3d 6e 2d 31 0d 00 64 | 36 e8 20 22 68 31 20 3d |%=n-1..d|6. "h1 =|
|000000b0| 20 22 20 68 31 20 3a f4 | 20 50 61 72 61 6d 65 74 | " h1 :.| Paramet|
|000000c0| 65 72 20 66 6f 72 20 4c | 69 6e 65 2d 53 70 72 65 |er for L|ine-Spre|
|000000d0| 61 64 2d 46 75 6e 63 74 | 69 6f 6e 0d 00 6e 3e 68 |ad-Funct|ion..n>h|
|000000e0| 32 3d 31 2d 68 31 20 20 | 20 20 20 20 20 20 20 20 |2=1-h1  |        |
|000000f0| 3a f4 20 62 79 20 74 68 | 65 20 73 75 6d 20 6f 66 |:. by th|e sum of|
|00000100| 20 66 75 6e 63 74 69 6f | 6e 20 6f 66 20 74 77 6f | functio|n of two|
|00000110| 20 67 61 75 73 73 69 6f | 6e 0d 00 78 31 e8 20 22 | gaussio|n..x1. "|
|00000120| 77 31 20 3d 20 22 20 77 | 31 20 3a f4 20 66 75 6e |w1 = " w|1 :. fun|
|00000130| 63 74 69 6f 6e 20 77 69 | 74 68 20 61 6d 70 6c 69 |ction wi|th ampli|
|00000140| 74 75 64 65 20 68 31 2c | 68 32 0d 00 82 10 e8 20 |tude h1,|h2..... |
|00000150| 22 77 32 20 3d 20 22 20 | 77 32 0d 00 8c 44 e8 20 |"w2 = " |w2...D. |
|00000160| 22 67 61 20 3d 20 22 20 | 67 61 20 3a f4 20 67 61 |"ga = " |ga :. ga|
|00000170| 20 3d 20 6e 6f 69 73 65 | 20 74 6f 20 73 69 67 6e | = noise| to sign|
|00000180| 61 6c 20 70 6f 77 65 72 | 20 64 65 6e 73 69 74 79 |al power| density|
|00000190| 20 72 61 74 69 6f 20 28 | 7e 30 2e 30 31 29 0d 00 | ratio (|~0.01)..|
|000001a0| 96 21 de 20 72 65 73 75 | 6c 74 25 20 31 37 36 2c |.!. resu|lt% 176,|
|000001b0| 6f 72 67 20 26 32 30 30 | 30 2c 76 61 6c 20 34 0d |org &200|0,val 4.|
|000001c0| 00 a0 4c de 20 78 72 65 | 61 6c 28 6e 2b 31 29 2c |..L. xre|al(n+1),|
|000001d0| 78 69 6d 61 67 28 6e 2b | 31 29 2c 78 73 28 6e 2b |ximag(n+|1),xs(n+|
|000001e0| 31 29 2c 78 28 6e 2b 31 | 29 2c 79 72 28 6e 30 25 |1),x(n+1|),yr(n0%|
|000001f0| 2c 6e 30 25 29 2c 79 69 | 28 6e 30 25 2c 6e 30 25 |,n0%),yi|(n0%,n0%|
|00000200| 29 2c 73 28 6e 29 2c 63 | 28 6e 29 0d 00 aa 04 0d |),s(n),c|(n).....|
|00000210| 00 b4 1f f4 20 2a 53 45 | 54 20 57 69 6c 64 24 50 |.... *SE|T Wild$P|
|00000220| 61 74 68 20 3c 4f 62 65 | 79 24 44 69 72 3e 0d 00 |ath <Obe|y$Dir>..|
|00000230| be 16 f4 20 2a 44 49 52 | 20 3c 57 69 6c 64 24 50 |... *DIR| <Wild$P|
|00000240| 61 74 68 3e 0d 00 c8 17 | ff 22 4c 4f 41 44 20 43 |ath>....|."LOAD C|
|00000250| 4f 44 45 20 22 2b c3 7e | 6f 72 67 0d 00 d2 04 0d |ODE "+.~|org.....|
|00000260| 00 dc 0e 6d 61 78 5f 69 | 6f 3d 6f 72 67 0d 00 e6 |...max_i|o=org...|
|00000270| 10 73 77 61 70 31 30 3d | 6f 72 67 2b 34 0d 00 f0 |.swap10=|org+4...|
|00000280| 0e 66 66 74 31 3d 6f 72 | 67 2b 38 0d 00 fa 0f 66 |.fft1=or|g+8....f|
|00000290| 66 74 32 3d 6f 72 67 2b | 31 32 0d 01 04 15 66 66 |ft2=org+|12....ff|
|000002a0| 74 5f 66 69 6c 74 65 72 | 3d 6f 72 67 2b 31 36 0d |t_filter|=org+16.|
|000002b0| 01 0e 10 66 66 74 31 69 | 3d 6f 72 67 2b 32 30 0d |...fft1i|=org+20.|
|000002c0| 01 18 10 66 66 74 32 69 | 3d 6f 72 67 2b 32 34 0d |...fft2i|=org+24.|
|000002d0| 01 22 04 0d 01 2c 36 e3 | 20 69 25 3d 30 20 b8 20 |."...,6.| i%=0 . |
|000002e0| 32 35 36 3a 70 3d af 2a | 69 25 2f 31 32 38 3a 63 |256:p=.*|i%/128:c|
|000002f0| 28 69 25 29 3d 9b 28 70 | 29 3a 73 28 69 25 29 3d |(i%)=.(p|):s(i%)=|
|00000300| b5 28 70 29 3a ed 20 69 | 25 0d 01 36 33 e8 22 45 |.(p):. i|%..63."E|
|00000310| 69 6e 6c 65 67 65 6e 20 | 64 65 72 20 44 61 74 65 |inlegen |der Date|
|00000320| 6e 64 69 73 6b 65 74 74 | 65 20 2d 20 49 6e 70 75 |ndiskett|e - Inpu|
|00000330| 74 2d 4e 61 6d 65 20 3a | 22 20 61 24 0d 01 40 09 |t-Name :|" a$..@.|
|00000340| 2a 4d 4f 55 2e 0d 01 4a | 0c 41 25 3d ad 28 61 24 |*MOU...J|.A%=.(a$|
|00000350| 29 0d 01 54 10 e3 20 6a | 25 3d 30 20 b8 20 6e 30 |)..T.. j|%=0 . n0|
|00000360| 25 0d 01 5e 10 e3 20 69 | 25 3d 30 20 b8 20 6e 30 |%..^.. i|%=0 . n0|
|00000370| 25 0d 01 68 12 79 72 28 | 6a 25 2c 69 25 29 3d 9a |%..h.yr(|j%,i%)=.|
|00000380| 23 41 25 0d 01 72 08 ed | 20 69 25 0d 01 7c 08 ed |#A%..r..| i%..|..|
|00000390| 20 6a 25 0d 01 86 08 d9 | 23 41 25 0d 01 90 34 e8 | j%.....|#A%...4.|
|000003a0| 22 45 69 6e 6c 65 67 65 | 6e 20 64 65 72 20 44 61 |"Einlege|n der Da|
|000003b0| 74 65 6e 64 69 73 6b 65 | 74 74 65 20 2d 20 4f 75 |tendiske|tte - Ou|
|000003c0| 74 70 75 74 2d 4e 61 6d | 65 20 3a 22 20 62 24 0d |tput-Nam|e :" b$.|
|000003d0| 01 9a 09 2a 4d 4f 55 2e | 0d 01 a4 04 0d 01 ae 0e |...*MOU.|........|
|000003e0| 41 25 3d 72 65 73 75 6c | 74 25 0d 01 b8 04 0d 01 |A%=resul|t%......|
|000003f0| c2 05 db 0d 01 cc 13 f1 | 20 22 46 46 74 20 6c 61 |........| "FFt la|
|00000400| 65 75 66 74 22 2c 0d 01 | d6 07 74 3d 91 0d 01 e0 |euft",..|..t=....|
|00000410| 04 0d 01 ea 2f d6 20 66 | 66 74 31 2c 41 25 2c 79 |..../. f|ft1,A%,y|
|00000420| 69 28 29 2c 79 72 28 29 | 2c 73 28 29 2c 63 28 29 |i(),yr()|,s(),c()|
|00000430| 2c 78 69 6d 61 67 28 29 | 2c 78 72 65 61 6c 28 29 |,ximag()|,xreal()|
|00000440| 0d 01 f4 12 f1 20 22 46 | 46 74 20 31 2e 22 2c 91 |..... "F|Ft 1.",.|
|00000450| 2d 74 0d 01 fe 2f d6 20 | 66 66 74 32 2c 41 25 2c |-t.../. |fft2,A%,|
|00000460| 79 69 28 29 2c 79 72 28 | 29 2c 73 28 29 2c 63 28 |yi(),yr(|),s(),c(|
|00000470| 29 2c 78 69 6d 61 67 28 | 29 2c 78 72 65 61 6c 28 |),ximag(|),xreal(|
|00000480| 29 0d 02 08 12 f1 20 22 | 46 46 74 20 32 2e 22 2c |)..... "|FFt 2.",|
|00000490| 91 2d 74 0d 02 12 04 0d | 02 1c 1a f4 20 41 6e 77 |.-t.....|.... Anw|
|000004a0| 65 6e 64 75 6e 67 20 64 | 65 72 20 46 69 6c 74 65 |endung d|er Filte|
|000004b0| 72 0d 02 26 04 0d 02 30 | 0b 6e 32 25 3d 6e 2f 32 |r..&...0|.n2%=n/2|
|000004c0| 0d 02 3a 04 0d 02 44 10 | e3 20 69 25 3d 31 20 b8 |..:...D.|. i%=1 .|
|000004d0| 20 6e 32 25 0d 02 4e 3a | 78 72 65 61 6c 28 69 25 | n2%..N:|xreal(i%|
|000004e0| 2b 6e 2f 32 29 3d a4 6c | 73 66 28 69 25 2d 31 29 |+n/2)=.l|sf(i%-1)|
|000004f0| 3a 78 72 65 61 6c 28 6e | 2f 32 2b 31 2d 69 25 29 |:xreal(n|/2+1-i%)|
|00000500| 3d a4 6c 73 66 28 69 25 | 29 3a ed 20 69 25 0d 02 |=.lsf(i%|):. i%..|
|00000510| 58 2c e3 20 69 25 3d 31 | 20 b8 20 6e 20 88 20 32 |X,. i%=1| . n . 2|
|00000520| 3a 78 72 65 61 6c 28 69 | 25 29 3d 2d 78 72 65 61 |:xreal(i|%)=-xrea|
|00000530| 6c 28 69 25 29 3a ed 20 | 69 25 0d 02 62 0d 78 69 |l(i%):. |i%..b.xi|
|00000540| 6d 61 67 28 29 3d 30 0d | 02 6c 2b d6 20 66 66 74 |mag()=0.|.l+. fft|
|00000550| 5f 66 69 6c 74 65 72 2c | 41 25 2c 73 28 29 2c 63 |_filter,|A%,s(),c|
|00000560| 28 29 2c 78 69 6d 61 67 | 28 29 2c 78 72 65 61 6c |(),ximag|(),xreal|
|00000570| 28 29 0d 02 76 09 6d 61 | 78 3d 30 0d 02 80 0e e3 |()..v.ma|x=0.....|
|00000580| 20 69 25 3d 31 20 b8 20 | 6e 0d 02 8a 44 78 73 28 | i%=1 . |n...Dxs(|
|00000590| 69 25 29 3d 28 b6 28 78 | 72 65 61 6c 28 69 25 29 |i%)=(.(x|real(i%)|
|000005a0| 5e 32 2b 78 69 6d 61 67 | 28 69 25 29 5e 32 29 29 |^2+ximag|(i%)^2))|
|000005b0| 2f 28 78 72 65 61 6c 28 | 69 25 29 5e 32 2b 78 69 |/(xreal(|i%)^2+xi|
|000005c0| 6d 61 67 28 69 25 29 5e | 32 2b 67 61 29 0d 02 94 |mag(i%)^|2+ga)...|
|000005d0| 23 e7 20 94 28 78 73 28 | 69 25 29 29 3e 6d 61 78 |#. .(xs(|i%))>max|
|000005e0| 20 8c 20 6d 61 78 3d 94 | 28 78 73 28 69 25 29 29 | . max=.|(xs(i%))|
|000005f0| 0d 02 9e 08 ed 20 69 25 | 0d 02 a8 25 e3 20 69 25 |..... i%|...%. i%|
|00000600| 3d 30 20 b8 20 6e 3a 78 | 73 28 69 25 29 3d 78 73 |=0 . n:x|s(i%)=xs|
|00000610| 28 69 25 29 2f 6d 61 78 | 3a ed 20 69 25 0d 02 b2 |(i%)/max|:. i%...|
|00000620| 04 0d 02 bc 22 f4 20 6c | 6f 77 20 66 72 65 71 75 |....". l|ow frequ|
|00000630| 65 6e 63 79 20 63 75 74 | 2d 6f 66 66 20 66 69 6c |ency cut|-off fil|
|00000640| 74 65 72 0d 02 c6 21 f4 | 20 20 20 20 20 70 6c 75 |ter...!.|     plu|
|00000650| 73 20 72 65 73 74 6f 72 | 61 74 69 6f 6e 20 66 69 |s restor|ation fi|
|00000660| 6c 74 65 72 0d 02 d0 04 | 0d 02 da 09 63 31 3d 2e |lter....|....c1=.|
|00000670| 36 0d 02 e4 0c 63 32 3d | 2e 30 30 37 31 0d 02 ee |6....c2=|.0071...|
|00000680| 14 70 3d 28 31 2d 63 31 | 29 2a 78 73 28 6e 32 25 |.p=(1-c1|)*xs(n2%|
|00000690| 29 0d 02 f8 1d 79 72 28 | 6e 32 25 2c 6e 32 25 29 |)....yr(|n2%,n2%)|
|000006a0| 3d 79 72 28 6e 32 25 2c | 6e 32 25 29 2a 70 0d 03 |=yr(n2%,|n2%)*p..|
|000006b0| 02 1d 79 69 28 6e 32 25 | 2c 6e 32 25 29 3d 79 69 |..yi(n2%|,n2%)=yi|
|000006c0| 28 6e 32 25 2c 6e 32 25 | 29 2a 70 0d 03 0c 12 e3 |(n2%,n2%|)*p.....|
|000006d0| 20 6a 25 3d 31 20 b8 20 | 6e 32 25 2d 31 0d 03 16 | j%=1 . |n2%-1...|
|000006e0| 2b 70 3d 28 31 2d 63 31 | 2a a1 28 2d 31 2a 28 6a |+p=(1-c1|*.(-1*(j|
|000006f0| 25 2f 6e 30 25 29 5e 32 | 2f 63 32 29 29 2a 78 73 |%/n0%)^2|/c2))*xs|
|00000700| 28 6e 32 25 2d 6a 25 29 | 0d 03 20 23 79 72 28 6e |(n2%-j%)|.. #yr(n|
|00000710| 32 25 2c 6e 32 25 2d 6a | 25 29 3d 79 72 28 6e 32 |2%,n2%-j|%)=yr(n2|
|00000720| 25 2c 6e 32 25 2d 6a 25 | 29 2a 70 0d 03 2a 23 79 |%,n2%-j%|)*p..*#y|
|00000730| 72 28 6e 32 25 2c 6e 32 | 25 2b 6a 25 29 3d 79 72 |r(n2%,n2|%+j%)=yr|
|00000740| 28 6e 32 25 2c 6e 32 25 | 2b 6a 25 29 2a 70 0d 03 |(n2%,n2%|+j%)*p..|
|00000750| 34 23 79 69 28 6e 32 25 | 2c 6e 32 25 2d 6a 25 29 |4#yi(n2%|,n2%-j%)|
|00000760| 3d 79 69 28 6e 32 25 2c | 6e 32 25 2d 6a 25 29 2a |=yi(n2%,|n2%-j%)*|
|00000770| 70 0d 03 3e 23 79 69 28 | 6e 32 25 2c 6e 32 25 2b |p..>#yi(|n2%,n2%+|
|00000780| 6a 25 29 3d 79 69 28 6e | 32 25 2c 6e 32 25 2b 6a |j%)=yi(n|2%,n2%+j|
|00000790| 25 29 2a 70 0d 03 48 23 | 79 72 28 6e 32 25 2d 6a |%)*p..H#|yr(n2%-j|
|000007a0| 25 2c 6e 32 25 29 3d 79 | 72 28 6e 32 25 2d 6a 25 |%,n2%)=y|r(n2%-j%|
|000007b0| 2c 6e 32 25 29 2a 70 0d | 03 52 23 79 72 28 6e 32 |,n2%)*p.|.R#yr(n2|
|000007c0| 25 2b 6a 25 2c 6e 32 25 | 29 3d 79 72 28 6e 32 25 |%+j%,n2%|)=yr(n2%|
|000007d0| 2b 6a 25 2c 6e 32 25 29 | 2a 70 0d 03 5c 23 79 69 |+j%,n2%)|*p..\#yi|
|000007e0| 28 6e 32 25 2d 6a 25 2c | 6e 32 25 29 3d 79 69 28 |(n2%-j%,|n2%)=yi(|
|000007f0| 6e 32 25 2d 6a 25 2c 6e | 32 25 29 2a 70 0d 03 66 |n2%-j%,n|2%)*p..f|
|00000800| 23 79 69 28 6e 32 25 2b | 6a 25 2c 6e 32 25 29 3d |#yi(n2%+|j%,n2%)=|
|00000810| 79 69 28 6e 32 25 2b 6a | 25 2c 6e 32 25 29 2a 70 |yi(n2%+j|%,n2%)*p|
|00000820| 0d 03 70 08 ed 20 6a 25 | 0d 03 7a 12 e3 20 6a 25 |..p.. j%|..z.. j%|
|00000830| 3d 31 20 b8 20 6e 32 25 | 2d 31 0d 03 84 13 e3 20 |=1 . n2%|-1..... |
|00000840| 70 72 25 3d 31 20 b8 20 | 6e 32 25 2d 31 0d 03 8e |pr%=1 . |n2%-1...|
|00000850| 2a 70 6a 3d b6 28 6a 25 | 5e 32 2b 70 72 25 5e 32 |*pj=.(j%|^2+pr%^2|
|00000860| 29 3a e7 20 70 6a 3e 6e | 32 25 2d 31 20 8c 20 70 |):. pj>n|2%-1 . p|
|00000870| 6a 3d 6e 32 25 2d 31 0d | 03 98 41 78 73 3d 28 78 |j=n2%-1.|..Axs=(x|
|00000880| 73 28 a8 28 6e 32 25 2d | 70 6a 29 29 2d 78 73 28 |s(.(n2%-|pj))-xs(|
|00000890| a8 28 6e 32 25 2b 31 2d | 70 6a 29 29 29 2a 28 70 |.(n2%+1-|pj)))*(p|
|000008a0| 6a 2d a8 28 70 6a 29 29 | 2b 78 73 28 a8 28 6e 32 |j-.(pj))|+xs(.(n2|
|000008b0| 25 2b 31 2d 70 6a 29 29 | 0d 03 a2 23 70 3d 28 31 |%+1-pj))|...#p=(1|
|000008c0| 2d 63 31 2a a1 28 2d 31 | 2a 28 70 6a 2f 6e 30 25 |-c1*.(-1|*(pj/n0%|
|000008d0| 29 5e 32 2f 63 32 29 29 | 2a 78 73 0d 03 ac 2b 79 |)^2/c2))|*xs...+y|
|000008e0| 72 28 6e 32 25 2d 70 72 | 25 2c 6e 32 25 2d 6a 25 |r(n2%-pr|%,n2%-j%|
|000008f0| 29 3d 79 72 28 6e 32 25 | 2d 70 72 25 2c 6e 32 25 |)=yr(n2%|-pr%,n2%|
|00000900| 2d 6a 25 29 2a 70 0d 03 | b6 2b 79 72 28 6e 32 25 |-j%)*p..|.+yr(n2%|
|00000910| 2b 70 72 25 2c 6e 32 25 | 2d 6a 25 29 3d 79 72 28 |+pr%,n2%|-j%)=yr(|
|00000920| 6e 32 25 2b 70 72 25 2c | 6e 32 25 2d 6a 25 29 2a |n2%+pr%,|n2%-j%)*|
|00000930| 70 0d 03 c0 2b 79 69 28 | 6e 32 25 2d 70 72 25 2c |p...+yi(|n2%-pr%,|
|00000940| 6e 32 25 2d 6a 25 29 3d | 79 69 28 6e 32 25 2d 70 |n2%-j%)=|yi(n2%-p|
|00000950| 72 25 2c 6e 32 25 2d 6a | 25 29 2a 70 0d 03 ca 2b |r%,n2%-j|%)*p...+|
|00000960| 79 69 28 6e 32 25 2b 70 | 72 25 2c 6e 32 25 2d 6a |yi(n2%+p|r%,n2%-j|
|00000970| 25 29 3d 79 69 28 6e 32 | 25 2b 70 72 25 2c 6e 32 |%)=yi(n2|%+pr%,n2|
|00000980| 25 2d 6a 25 29 2a 70 0d | 03 d4 2b 79 72 28 6e 32 |%-j%)*p.|..+yr(n2|
|00000990| 25 2d 70 72 25 2c 6e 32 | 25 2b 6a 25 29 3d 79 72 |%-pr%,n2|%+j%)=yr|
|000009a0| 28 6e 32 25 2d 70 72 25 | 2c 6e 32 25 2b 6a 25 29 |(n2%-pr%|,n2%+j%)|
|000009b0| 2a 70 0d 03 de 2b 79 72 | 28 6e 32 25 2b 70 72 25 |*p...+yr|(n2%+pr%|
|000009c0| 2c 6e 32 25 2b 6a 25 29 | 3d 79 72 28 6e 32 25 2b |,n2%+j%)|=yr(n2%+|
|000009d0| 70 72 25 2c 6e 32 25 2b | 6a 25 29 2a 70 0d 03 e8 |pr%,n2%+|j%)*p...|
|000009e0| 2b 79 69 28 6e 32 25 2d | 70 72 25 2c 6e 32 25 2b |+yi(n2%-|pr%,n2%+|
|000009f0| 6a 25 29 3d 79 69 28 6e | 32 25 2d 70 72 25 2c 6e |j%)=yi(n|2%-pr%,n|
|00000a00| 32 25 2b 6a 25 29 2a 70 | 0d 03 f2 2b 79 69 28 6e |2%+j%)*p|...+yi(n|
|00000a10| 32 25 2b 70 72 25 2c 6e | 32 25 2b 6a 25 29 3d 79 |2%+pr%,n|2%+j%)=y|
|00000a20| 69 28 6e 32 25 2b 70 72 | 25 2c 6e 32 25 2b 6a 25 |i(n2%+pr|%,n2%+j%|
|00000a30| 29 2a 70 0d 03 fc 09 ed | 20 70 72 25 0d 04 06 08 |)*p.....| pr%....|
|00000a40| ed 20 6a 25 0d 04 10 15 | f1 22 46 46 54 2d 46 69 |. j%....|."FFT-Fi|
|00000a50| 6c 74 65 72 22 2c 91 2d | 74 0d 04 1a 04 0d 04 24 |lter",.-|t......$|
|00000a60| 30 d6 20 66 66 74 31 69 | 2c 41 25 2c 79 69 28 29 |0. fft1i|,A%,yi()|
|00000a70| 2c 79 72 28 29 2c 73 28 | 29 2c 63 28 29 2c 78 69 |,yr(),s(|),c(),xi|
|00000a80| 6d 61 67 28 29 2c 78 72 | 65 61 6c 28 29 0d 04 2e |mag(),xr|eal()...|
|00000a90| 13 f1 20 22 46 46 74 20 | 31 69 2e 22 2c 91 2d 74 |.. "FFt |1i.",.-t|
|00000aa0| 0d 04 38 30 d6 20 66 66 | 74 32 69 2c 41 25 2c 79 |..80. ff|t2i,A%,y|
|00000ab0| 69 28 29 2c 79 72 28 29 | 2c 73 28 29 2c 63 28 29 |i(),yr()|,s(),c()|
|00000ac0| 2c 78 69 6d 61 67 28 29 | 2c 78 72 65 61 6c 28 29 |,ximag()|,xreal()|
|00000ad0| 0d 04 42 13 f1 20 22 46 | 46 74 20 32 69 2e 22 2c |..B.. "F|Ft 2i.",|
|00000ae0| 91 2d 74 0d 04 4c 04 0d | 04 56 2b e3 20 69 25 3d |.-t..L..|.V+. i%=|
|00000af0| 30 20 b8 20 6e 30 25 3a | 79 72 28 30 2c 69 25 29 |0 . n0%:|yr(0,i%)|
|00000b00| 3d 30 3a 79 72 28 31 2c | 69 25 29 3d 30 3a ed 20 |=0:yr(1,|i%)=0:. |
|00000b10| 69 25 0d 04 60 0a 7c 76 | 61 6c 3d 30 0d 04 6a 35 |i%..`.|v|al=0..j5|
|00000b20| 72 65 73 75 6c 74 25 21 | 31 36 38 3d 76 61 6c 21 |result%!|168=val!|
|00000b30| 30 3a 72 65 73 75 6c 74 | 25 21 31 37 32 3d 30 3a |0:result|%!172=0:|
|00000b40| 72 65 73 75 6c 74 25 3f | 31 37 32 3d 76 61 6c 3f |result%?|172=val?|
|00000b50| 34 0d 04 74 14 d6 20 6d | 61 78 5f 69 6f 2c 41 25 |4..t.. m|ax_io,A%|
|00000b60| 2c 79 72 28 29 0d 04 7e | 31 76 61 6c 21 30 3d 72 |,yr()..~|1val!0=r|
|00000b70| 65 73 75 6c 74 25 21 31 | 36 38 3a 76 61 6c 3f 34 |esult%!1|68:val?4|
|00000b80| 3d 72 65 73 75 6c 74 25 | 3f 31 37 32 3a 6d 61 78 |=result%|?172:max|
|00000b90| 6f 3d 7c 76 61 6c 0d 04 | 88 04 0d 04 92 13 6d 61 |o=|val..|......ma|
|00000ba0| 78 72 3d 28 6e 2d 31 29 | 2f 6d 61 78 6f 0d 04 9c |xr=(n-1)|/maxo...|
|00000bb0| 14 d6 20 73 77 61 70 31 | 30 2c 41 25 2c 79 72 28 |.. swap1|0,A%,yr(|
|00000bc0| 29 0d 04 a6 04 0d 04 b0 | 09 74 3d 91 2d 74 0d 04 |).......|.t=.-t..|
|00000bd0| ba 0c 41 25 3d ae 28 62 | 24 29 0d 04 c4 11 e3 20 |..A%=.(b|$)..... |
|00000be0| 70 72 25 3d 30 20 b8 20 | 6e 30 25 0d 04 ce 10 e3 |pr%=0 . |n0%.....|
|00000bf0| 20 69 25 3d 30 20 b8 20 | 6e 30 25 0d 04 d8 14 e7 | i%=0 . |n0%.....|
|00000c00| 20 79 72 28 70 72 25 2c | 69 25 29 3e 30 20 8c 0d | yr(pr%,|i%)>0 ..|
|00000c10| 04 e2 1e d5 23 41 25 2c | a8 28 79 72 28 70 72 25 |....#A%,|.(yr(pr%|
|00000c20| 2c 69 25 29 2a 6d 61 78 | 72 2b 2e 35 29 0d 04 ec |,i%)*max|r+.5)...|
|00000c30| 05 cc 0d 04 f6 0a d5 23 | 41 25 2c 30 0d 05 00 05 |.......#|A%,0....|
|00000c40| cd 0d 05 0a 08 ed 20 69 | 25 0d 05 14 09 ed 20 70 |...... i|%..... p|
|00000c50| 72 25 0d 05 1e 08 d9 23 | 41 25 0d 05 28 0a f2 70 |r%.....#|A%..(..p|
|00000c60| 72 69 6e 74 0d 05 32 05 | e0 0d 05 3c 04 0d 05 46 |rint..2.|...<...F|
|00000c70| 0b dd f2 70 72 69 6e 74 | 0d 05 50 0f db 3a da 3a |...print|..P..:.:|
|00000c80| f1 20 74 2f 31 30 30 0d | 05 5a 1e e8 22 58 2c 59 |. t/100.|.Z.."X,Y|
|00000c90| 20 2d 20 45 20 3a 20 4e | 72 2e 20 3a 20 22 20 61 | - E : N|r. : " a|
|00000ca0| 24 2c 6e 72 25 0d 05 64 | 13 e7 20 61 24 3d 22 45 |$,nr%..d|.. a$="E|
|00000cb0| 22 20 8c 20 8d 64 7a 45 | 0d 05 6e 0e e7 20 61 24 |" . .dzE|..n.. a$|
|00000cc0| 3d 22 59 22 20 8c 0d 05 | 78 10 e3 20 69 25 3d 30 |="Y" ...|x.. i%=0|
|00000cd0| 20 b8 20 6e 30 25 0d 05 | 82 1f 78 72 65 61 6c 28 | . n0%..|..xreal(|
|00000ce0| 69 25 2b 31 29 3d 79 72 | 28 69 25 2c 6e 72 25 29 |i%+1)=yr|(i%,nr%)|
|00000cf0| 2a 6d 61 78 72 0d 05 8c | 1f 78 69 6d 61 67 28 69 |*maxr...|.ximag(i|
|00000d00| 25 2b 31 29 3d 79 69 28 | 69 25 2c 6e 72 25 29 2a |%+1)=yi(|i%,nr%)*|
|00000d10| 6d 61 78 72 0d 05 96 08 | ed 20 69 25 0d 05 a0 05 |maxr....|. i%....|
|00000d20| cc 0d 05 aa 10 e3 20 69 | 25 3d 30 20 b8 20 6e 30 |...... i|%=0 . n0|
|00000d30| 25 0d 05 b4 1f 78 72 65 | 61 6c 28 69 25 2b 31 29 |%....xre|al(i%+1)|
|00000d40| 3d 79 72 28 6e 72 25 2c | 69 25 29 2a 6d 61 78 72 |=yr(nr%,|i%)*maxr|
|00000d50| 0d 05 be 1f 78 69 6d 61 | 67 28 69 25 2b 31 29 3d |....xima|g(i%+1)=|
|00000d60| 79 69 28 6e 72 25 2c 69 | 25 29 2a 6d 61 78 72 0d |yi(nr%,i|%)*maxr.|
|00000d70| 05 c8 08 ed 20 69 25 0d | 05 d2 05 cd 0d 05 dc 12 |.... i%.|........|
|00000d80| f2 73 68 6f 77 28 78 72 | 65 61 6c 28 29 29 0d 05 |.show(xr|eal())..|
|00000d90| e6 12 f2 73 68 6f 77 28 | 78 69 6d 61 67 28 29 29 |...show(|ximag())|
|00000da0| 0d 05 f0 0a e5 20 8d 44 | 50 45 0d 05 fa 05 e1 0d |..... .D|PE......|
|00000db0| 06 04 04 0d 06 0e 0f dd | f2 73 68 6f 77 28 78 28 |........|.show(x(|
|00000dc0| 29 29 0d 06 18 07 db 3a | da 0d 06 22 18 e8 22 50 |)).....:|...".."P|
|00000dd0| 72 69 6e 74 65 72 20 28 | 4a 2f 4e 29 20 22 2c 61 |rinter (|J/N) ",a|
|00000de0| 24 0d 06 2c 12 e7 20 61 | 24 3d 22 4a 22 20 8c 20 |$..,.. a|$="J" . |
|00000df0| ef 20 32 0d 06 36 0e e3 | 20 69 25 3d 31 20 b8 20 |. 2..6..| i%=1 . |
|00000e00| 6e 0d 06 40 20 e7 20 28 | 69 25 20 83 20 36 29 3d |n..@ . (|i% . 6)=|
|00000e10| 31 20 8c 20 40 25 3d 26 | 33 30 34 3a f1 20 69 25 |1 . @%=&|304:. i%|
|00000e20| 3b 0d 06 4a 0b 40 25 3d | 26 36 30 43 0d 06 54 0c |;..J.@%=|&60C..T.|
|00000e30| f1 20 78 28 69 25 29 3b | 0d 06 5e 14 e7 20 28 69 |. x(i%);|..^.. (i|
|00000e40| 25 20 83 20 36 29 3d 30 | 20 8c 20 f1 0d 06 68 08 |% . 6)=0| . ...h.|
|00000e50| ed 20 69 25 0d 06 72 05 | f1 0d 06 7c 07 ef 20 33 |. i%..r.|...|.. 3|
|00000e60| 0d 06 86 08 e8 20 61 24 | 0d 06 90 07 db 3a da 0d |..... a$|.....:..|
|00000e70| 06 9a 09 6d 61 78 3d 30 | 0d 06 a4 0e e3 20 69 25 |...max=0|..... i%|
|00000e80| 3d 31 20 b8 20 6e 0d 06 | ae 21 e7 20 94 28 78 28 |=1 . n..|.!. .(x(|
|00000e90| 69 25 29 29 3e 6d 61 78 | 20 8c 20 6d 61 78 3d 94 |i%))>max| . max=.|
|00000ea0| 28 78 28 69 25 29 29 0d | 06 b8 08 ed 20 69 25 0d |(x(i%)).|.... i%.|
|00000eb0| 06 c2 12 f1 20 22 6d 61 | 78 20 3d 20 22 3b 6d 61 |.... "ma|x = ";ma|
|00000ec0| 78 0d 06 cc 13 e7 20 6d | 61 78 3d 30 20 8c 20 6d |x..... m|ax=0 . m|
|00000ed0| 61 78 3d 31 0d 06 d6 07 | e6 20 33 0d 06 e0 0c c8 |ax=1....|. 3.....|
|00000ee0| 91 20 30 2c 35 31 31 0d | 06 ea 10 86 20 30 2c 30 |. 0,511.|.... 0,0|
|00000ef0| 2c 31 30 32 34 2c 30 0d | 06 f4 07 e6 20 31 0d 06 |,1024,0.|.... 1..|
|00000f00| fe 14 ec 20 30 2c 78 28 | 31 29 2a 35 31 32 2f 6d |... 0,x(|1)*512/m|
|00000f10| 61 78 0d 07 08 0f e3 20 | 69 3d 31 20 b8 20 32 35 |ax..... |i=1 . 25|
|00000f20| 36 0d 07 12 1c f0 20 35 | 2c 28 69 2d 31 29 2a 34 |6..... 5|,(i-1)*4|
|00000f30| 2c 78 28 69 29 2a 35 31 | 32 2f 6d 61 78 0d 07 1c |,x(i)*51|2/max...|
|00000f40| 07 ed 20 69 0d 07 26 08 | e8 20 61 24 0d 07 30 05 |.. i..&.|. a$..0.|
|00000f50| e1 0d 07 3a 04 0d 07 44 | 0c dd a4 6c 73 66 28 6d |...:...D|...lsf(m|
|00000f60| 29 0d 07 4e 1f e7 20 68 | 32 3d 30 20 80 20 6d 3e |)..N.. h|2=0 . m>|
|00000f70| 31 2e 35 2a 77 31 20 8c | 20 6d 3d 6d 2a 31 2e 35 |1.5*w1 .| m=m*1.5|
|00000f80| 0d 07 58 2b 6c 73 66 3d | 68 31 2a a1 28 2d 32 2a |..X+lsf=|h1*.(-2*|
|00000f90| 28 6d 2f 77 31 29 5e 32 | 29 2b 68 32 2a a1 28 2d |(m/w1)^2|)+h2*.(-|
|00000fa0| 32 2a 28 6d 2f 77 32 29 | 5e 32 29 0d 07 62 08 3d |2*(m/w2)|^2)..b.=|
|00000fb0| 6c 73 66 0d 07 6c 04 0d | 07 76 05 e0 0d ff       |lsf..l..|.v....  |
+--------+-------------------------+-------------------------+--------+--------+